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Abstract 

Piwi proteins are part of a superfamily of Argonaute proteins, wliicli are one of the core components of tine RNA silencing 
patliway in many eul<aryotes. Piwi proteins are tliought to repress the transposon expression both transcriptionally and 
post-transcriptionally. Recently, Drosophila melanogaster Piwi was recently reported to associate with chromatin and to 
interact directly with the Heterochromatin Protein 1 (HPla). However, similar interactions have not been reported in other 
higher eukaryotes. Here we show that silkworm Piwi proteins interact with HPIs in the nucleus. The silkworm, Bombyx mori, 
has two Piwi proteins, Ago3 and Siwi, and two typical HP1 proteins, HPla and HPl b. We found that HPla plays an important 
role in the interaction between Ago3/Siwi and HPl b in the ovary-derived BmN4 cell line. We also found that Ago3/Siwi 
regulates the transcription in an HPl -dependent manner. These results suggest that silkworm Piwi proteins function as a 
chromatin regulator in collaboration with HPla and HPIb. 
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Introduction 

RNA silencing triggered by 20-35-nucleotide (nt) small RNA is 
involved in the regulation of gene expression both transcriptionally 
and post-transcriptionaUy [1-3]. The RNA silencing system in the 
fission yeast, Schizosaccharomyces pombe, is particularly well charac- 
terized [2,3]. In S. pombe, which does not have Piwi-interacting 
RNAs (piRNAs), small interfering RNAs (siRNAs) and Argonaute 
1 (Agol) recognize nascent transcripts and assemble centromeric 
heterochromatin by recruiting the Heterochromatin Protein 1 
(HPl) homolog Swi6 [2,4]. In higher eukaryotes, multiple 
Argonaute proteins have been found and divided into two 
subfamilies [5]. The Piwi subfamily proteins, part of the 
superfamily of Argonaute proteins, are also thought to have an 
important role in the regnlation of heterochromatin and chroma- 
tin maintenance [1,6-8]. However, the silencing mechanism 
mediated by Piwi proteins has not been well understood. 

Drosophila melanogaster Piwi (dPiwi) was recently reported to 
associate with chromatin and to interact directly with D. 
melanogaster HPla (dHPla). Thus, it is possible that dPiwi silences 
transcription by heterochromatin assembly analogous to the case 
in S. pombe [7-9]. HPl has been generally known to be a non- 
histone chromosomal protein that plays diverse and critical roles in 
chromatin structure, transcription, DNA replication, chromosome 
segregation, and genomic stability [3]. The interaction between 
dPiwi and dHPla is mediated by a chromoshadow domain in 
dHP 1 a and a PxVxL-type motif in dPiwi [9] . In addition to the 
direct interaction between dPiwi and dHPla, dPiwi has been 
reported to guide epigenetic factors to their target sites and to be 



required for the transcriptional silencing of transposons in 
Drosophila [7,8]. 

The silkworm, Bombyx mori, has two Piwi genes encoding Ago3 
and Siwi, and two typical HPl proteins, HPla and HP lb, in the 
genome [10]. Similarly to Drosophila, silkworm Ago3 and Siwi play 
an important role, at least in part, in repressing transposons 
[11,12]. We previously reported that silkworm HPla forms a 
homo- and heterodimer with HP lb and interacts more strongly 
with Su(var)3-9 than does HP lb, whereas HP lb has higher 
repression activity than HPla [10]. Although direct interaction 
between Piwi and HPla proteins has been reported in Drosophila, 
similar interactions have not been reported in other higher 
eukaryotes [9], to our knowledge. Here, we report the functional 
relationship between Piwi proteins and HPl in the silkworm 
ovary-derived cultured cell line BmN4. The silkworm Piwi 
proteins were found to interact with HP lb through the mediation 
of HPla, and Ago3/Siwi was shown to be involved in 
transcriptional regulation in an HPl -dependent manner. These 
findings suggest that Ago3/Siwi may function as the chromatin 
regulator in silkworms. 

Results and Discussion 

Silkworm Piwi Proteins Interact with HPl in Nucleus 

To test the direct interaction between dPiwi and dHPla 
proteins in the silkworm ovary-derived cultured cell line BmN4, 
we employed BiFC (bimolecular fluorescence complementation; 
[13]) and analyzed the proteins' localization. The signal of dPiwi- 
dHPla interaction was observed primarily in the nucleus of the 
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BmN4 cells (Figure SI A), whereas the EGFP-fused dPiwi was 
diffused in BmN4 cells as opposed to Drosophila cells ([14]; Figure 
SIB). Using this BiFC assay system, we examined the interactions 
of two silkworm Piwi-subfamily proteins, Ago3 and Siwi, with two 
canonical silkworm HPl proteins, HP la and -b, in BmN4 cells. As 
shown in Figure lA, both Ago3 and Siwi interacted with both 
HPla and HPlb primarily in the nucleus of the BmN4 cells. 
Interestingly, the interaction signals were also observed in tiny 
cytoplasmic granule-like structures in a certain number of BmN4 
cells. Previously, we found similar Ago3 or Siwi granules in BmN4 
cells [1 1], but it is not yet clear whether the granule-like structures 
observed in the present study are the same as those observed 
previously. 

Unfortunately, a co-immunoprecipitation assay could not verify 
these interactions in BmN4 cells (data not shown), suggesting that 
silkworm Piwi interacts loosely or transiently with HPl proteins. 
Thus, we carried out an insect two-hybrid (I2H) analysis, a method 
with high sensitivity for the detection of this weak protein-protein 
interaction [15]. Similarly to the results obtained with BiFC, the 
co-expression of DBD-Ago3 or -Siwi (as bait) with either of the two 
AD-HP 1 proteins (as prey) increased the reporter activity more 
than 10-fold in BmN4 cells (Figure IB and C). The I2H system 
uses the GAL4 DNA-binding domain (GAL4 DBD) with a typical 
nuclear localization signal (NLS); thus, there remains the 
possibility that the DBD-mediated nuclear import of silkworm 
Piwi proteins may lead to false-positi\'e results for the interaction 
between silkworm Piwi and HPl proteins. 

To exclude this possibility, we constructed another I2H system 
in which TetR without typical NLS and tetO sequences was used 
instead of the GAL4 DBD and UAS sequences, respectively. In 
agrccmc-nt with the r(tsults obtained with the GAL4-based I2H, 
the reporter activity of the TetR-based I2H analysis indicated the 
interaction between silkworm HPl and Piwi-subfamUy proteins 
(Figures SIC and D). Taken together, these results clearly indicate 
that silkworm Piwi proteins interact with two canonical HPl 
proteins in the nucleus of BmN4 cultured cells. 

In Drosophila, dimeric dHP 1 a binds to a PxVxL-type motif in the 
N-terminal domain of dPiwi, which is found only in the Drosophila 
Piwi protein, through its C-terminal chromoshadow domain [9]. 
To investigate which region of the Piwi protein interacts with 
HPls, i.e., the N- or C-terminal domain, we performed I2H 
analyses between split silkworm Piwi and HPl proteins (Figure 
S2). In all silkworm Piwi proteins and dPiwi undergoing the I2H 
assay with HPl proteins, the N-terminal domain showed stronger 
luciferase activities than the C-terminal domain (Figure S3). Thus, 
similar to the case in Drosophila, these results suggest that the N- 
terminal domains of Piwi have a major contribution to the 
interaction with HPla/b in silkworm BmN4 cells. 

The silkworm genome sequence database, however, revealed 
that Ago3 and Siwi do not contain PxVxL-type motifs in their 
amino acid sequence, suggesting that the silkworm Piwi-subfamily 
proteins interact with HPls at regions distinct from the PxVxL- 
type motif, unlike other canonical HPl interactors [16]. 
Additionally, the interaction of the N-terminal domain of Piwi 
proteins with HPls was stronger than that of fiiU-length Piwi with 
HPls, even in the case of dPiwi-dHPla interaction (Figure S3). 
This result may indicate that the potential high affinities of N- 
terminal domains for I IP Is are modulated by other regions of the 
Piwi proteins. These results thus suggest the possibility that Piwi 
proteins of other species lacking typical HPl interaction motifs can 
interact with HPl in the same way as silkworm Piwi proteins. 



HPla Mediates the Interaction between Silkworm Piwi 
Proteins and HPlb 

Since silkworm HPla and HPlb have been demonstrated to 
form a heterodimeric complex [10], we analyzed the mutual 
dependency of HPls in the interaction with Ago3 and Siwi. Using 
BmN4-SIDl cells, a cell line highly sensitive to soaking RNAi [17], 
we performed I2H analyses under knockdown conditions of each 
HPl protein. Semi-quantitative RT-PCR analyses clearly demon- 
strated the specific knockdown of HPl mRNAs by RNAi 
(Figure 2A). Compared to the dsGFP-treated cells used as a 
control, the interaction of HPlb with Ago3 or Siwi was decreased 
in the absence of HPla (Figure 2C), suggesting that the interaction 
between Ago3/Siwi and HPlb is mediated by HPla (Figure 2D). 

In contrast, the interaction of HPla with Ago3 or Siwi detected 
by I2H was increased in BmN4 cells with down-regulated HPlb 
expression (Figure 2B). As illustrated in Figure 2E, the increased 
HPla homodimer formation in the absence of HPlb would 
enhance its interaction with Ago3 and Siwi, and the exclusion of 
the HPla/b heterodimer would remove the negative transcrip- 
tional activity of HPlb. 

We previously reported that the silkworm Ago3 and Siwi co- 
localize in nuage-like cytoplasmic granules of BmN4 cells [11]. In 
Drosophila, two piRNA-related cytoplasmic granules, nuages in 
nurse cells and Yb bodies in ovarian somatic cells (OSCs), have 
been reported [14,18]. In nurse cells, the localization of piRNA- 
related proteins to nuages depends on Vasa, a marker for germ 
cells [18]. In contrast, Armitage and Yb, the Yb body components, 
are required for the entry of dPiwi into the OSC nucleus [14]. We 
thus examined the effect of depletion of these piRNA-related genes 
on the interaction between Ago3/Siwi and HPla/b. As nuage- 
related genes, we used SPJ^-E and HID in addition to VLG (V asa- 
like gene; [19]). There were weak but not significant effects in the 
interaction between Ago3/Siwi and HPla/b under knockdown 
conditions (Figures S4 and S5), suggesting that Ago3 and Siwi 
interact with HPla independently of other piRNA-related genes. 
Only the Siwi-HPla interaction was increased significantiy in 
VLG knockdown cells (Figure S5C). 

In agreement with our previous findings [19], the present BiFC 
analysis detected the interaction of Ago3-Siwi in nuage-like 
granules, which are marked by the localization of DsRed-fused 
VTG (Figure 3). Interestingly, the Ago3-Siwi signal increased and 
became clearer under the HPla or -b depletion conditions, 
suggesting that HPl proteins stimulate the nuclear import of Ago3 
or Siwi from nuage-hke granules. Taken together, the past and 
present findings indicate that silkworm HPls may compete or 
coordinate with VLG to determine the subcellular localization of 
Piwi proteins predominantiy through Siwi-HPla interaction. 

Roles of Silkworm Piwi Proteins in Transcription 

On the basis of the existence of a Piwi-HPl interaction, Ago3/ 
Siwi might contribute to transcriptional gene silencing (TGS) 
through the induction of heterochromatin formation as do S. pombe 
Agol and D. melanogaster Piwi [1,2,4,6-9,20]. To assess whether 
Ago3/Siwi could function as a transcriptional repressor, we 
performed a transcription repression assay using the GAL4-UAS- 
based assay system previously reported [10]- We found that GAL4 
DBD-fused Piwi proteins could be forcibly recruited to UAS 
elements upstream of the reporter luciferase gene without 
guidance by piRNA [10]. As previously reported, HPl proteins 
repressed the transcription, and the repression induced by DBD- 
HPlb was more efficient than that by DBD-HPla. Compared to 
HPl proteins, DBD-Ago3 and -Siwi have only modest repression 
activity (Figure 4A). 
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Bait: DBD-Ago3 Bait: DBD-Siwi 



Figure 1. Silkworm Piwi proteins interacted with HP1a and -b in cultured silkworm cells. (A) The BiFC analysis of the interaction between 
B. mori Piwi and HP1 proteins in the silkworm cultured cell line BmN4. Venus and Hoechst33342 fluorescences indicate their localization and nucleus, 
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respectively. Scale bar: 10 nm. (B, C) The insect two-hybrid (I2H) assay of the interaction between silkworm Piwi and HP1 proteins. BmN4 cells were 
transfected with 4xUAS-TATA-Luc reporter construct, along with expression vectors for GAL4 DNA-binding domain (DBD) and p65 activation domain 
(AD) fused to Ago3, Siwi, HPIa or HPIb as indicated. Empty cloning vector pUC18 was introduced as the control prey instead of pAD-HPIa or HPlb. 
Luciferase activities were measured at 72 h post-transfection. Values are relative luciferase activities in the pUC18-introduced control. Error 
bars = standard deviation (SD). The SDs and P-values are based on n = 3, and were determined by the t-test in comparison with the luciferase activities 
in the pUC18-introduced controls. *P<0.05, **P<0.01, ***P<0.001. 
doi:1 0.1 371/journal.pone.009231 3.g001 
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Figure 2. HPIa is required for the interaction between silkworm Piwi proteins and HPIb. (A) Semi-quantitative RT-PCR analyses were 
performed using the cDNA library from AgoB-, Slwl-, HPIa- or HP1 b-gene-knockdown BmN4-SID1 cells. RNAI-lnduced BmN4-SID1 cells showed a 
reduction of target gene expression level 72 h after dsRNA introduction. (B, C) The I2H assay of the interaction between silkworm Plwl and HP1 
proteins In each HP1 -depleted BmN4-SID1 cell. At 72 h after the Introduction of dsHPIa, dsHPlb or dsGFP (control), BmN4-SID cells were transfected 
with 4xUAS-TATA-Luc reporter construct, along with expression vectors for GAL4 DNA-bindIng domain (DBD) and p65 activation domain (AD) fused 
to Ago3, Slwl, HPIa or -b as Indicated. Luciferase activities were measured at 72 h post-transfectlon. Luciferase activities were normalized to that of 
dsGFP-lntroduced controls (n = 3). Error bars = SD. The SDs and P-values (t-test, *P<0.05) are based on n = 3, and are based on comparisons with 
luciferase activities In dsGFP-lntroduced controls. (D, E) A model for the interaction between silkworm Piwi and HP1 proteins in each HP1 -gene- 
knockdown BmN4-SID1 cell. 
doi:1 0.1 371/journal.pone.009231 3.g002 
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Figure 3. The Ago3-Siwi signal in nuage-like granules increased under the HPIa or -b knocked-down conditions in the silkworm 
cultured cells. The BiFC analyses for the interaction between Ago3 and Siwi proteins in each HP1 -gene-knockdown BmN4-SID1 cell. The BiFC 
analyses were performed at 72 h after the introduction of dsHPIa, dsHPb or dsLuc (control). Venus, DsRed and 4',6-diamidino-2-phenylindole (DAPI) 
fluorescences indicate Ago3-Siwi, VLG (as a marker for nuage) and the nucleus, respectively. Scale bar: 10 |am. 
doi:1 0.1 371 /journal.pone.009231 3.g003 



We next attempted to examine the repression activity of 
silkworm Piwi proteins in BmN4-SIDl cells under conditions in 
wliich HPla, HPlb, Ago3 or Siwi were knocked down (Figure 4). 
The transcriptional repression mediated by DBD-Ago3 and -Siwi 
in the absence of HPla was slighdy decreased compared to that of 
the control (Figures 4A and B), suggesting that Ago3 and Siwi by 
themselves have no function in the transcriptional repression. 
Interestingly, DBD-Ago3 or -Siwi enhanced the transcription in 
the absence of HPlb (Figure 4C), implying that Ago3/Siwi has 
another interactor which regulates transcription positively and 
competes with HPlb in binding to Ago3/Siwi. 

Indeed, it has been reported that small RNA machinery 
interacts with RNA polymerase II, the second-largest RNA 
subunit, and affects heterochromatic silencing in D. melanogaster 
and S. pombe [21,22]. These fmdings suggest that HPla mediates 
the interactions of Ago3/Siwi with both HPlb and a positive 
regulatory factor. Unlike the knockdown of HPla or HPlb, the 
transcription-repression activities of DBD-HPla and -HPlb were 
not aflFected by the presence or absence of silkworm Piwi genes 
(Figures 4D and E). In summary, we suspect that HPla functions 
as a "hub" for the epigenetic chromatin regulation. 

The chromatin-loading mechanism of sUkworm Piwi proteins at 
their target loci has not yet been elucidated, but is thought to be 



mediated by piRNA as in the S. pombe system [4]. The HPl 
proteins would be recruited to a target site in a sequence-specific 
manner guided by the Piwi-piRNA complex without histone 
modification, such as the typical heterochromatin maker H3K9 
methylation [3]. At the target site, Piwi-HPl complex functions as 
a rapid transcriptional repressor depending on the transcriptional 
repression activity of HPl proteins before histone modifications 
against selfish genetic elements are established. Namely, the 
silkworm Piwi complex with an HPla homodimer is a moderate 
transcriptional modulator and can recruit transcriptional activa- 
tors or epigenetic factors, such as Su(var)3-9 [23]. At an early 
stage of the fist-time invasion of selfish genetic elements, the 
invasion site of which might be epigenetically naked and cannot be 
recognized by HPl proteins, a certain amount of Piwi complex 
assembles with HP 1 a-HP 1 b heterodimers at the transposed selfish 
genetic elements guided by piRNA, represses their transcription 
and then forms a broad heterochromatic region by the ability of 
HPl proteins to spread. This proposed mechanism might be the 
first step of de novo heterochromatin formation on selfish genetic 
elements. 
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Control (dsGFP) HP1 a knockdown HPIb knockdown 





Figure 4. Silkworm Piwi proteins repressed transcription in an IHPI -dependent manner. (A-F) Transcriptional repression of luciferase 2P 
(luc2P) reporter gene under the control of the hsp promoter by silkworm Piwi and HP1 proteins in BmN4-SID1 cells with dsRNA against Ago3, Siwi, 
HPl a or HP1 b. The hsp promoter-based luc2P reporter contains five synthetic GAL4-UAS sites and was used in all transfections. BmN4-SID1 cells were 
transfected with a Luc2P reporter under the control of the hsp promoter, upstream of which five UAS were located for the hsp promoter. Cells were 
co-transfected with an expression vector for GAL4 DNA-binding domain (DBD)-fused Ago3, Siwi, HPl a or HPIb proteins at 72 h after dsRNA 
introduction. Luciferase activities were measured at 72 h post-transfection. Luciferase activities are represented as relative values based on DBD 
(empty plasmid introduced cells) as a standard. Error bars = SD. The SDs and P-values (determined by the t-test, *P<0.1, **P<0.05, ***P<0.01, which 
was checked by a comparison with the luciferase activities in the DBD-introduced controls with DBD-HPIa, DBD-HPIb, DBD-Ago3 or DBD-Siwi- 
introduced cells) are based on n = 3. 
dor:l 0.1 371/journal.pone.009231 3.g004 



Materials and Methods 

Construction of Plasmids 

The entry clones for dHPla and dPiwi were constructed on 
pENTRl 1 from Invitrogen. To obtain the entr)' clone for dHPla, 
the total RNA was isolated from the Dwsophila melanogaster 
Schneider 2 cultured cell. This total RNA was reverse-transcribed 
with oUgo (dT) primer and amplified with the primer set: 5'- 
GGGCC ATGGG CAAGA AAATC GACAA CCCTG-3' and 
5'-GGGCT CGAGT TAATC TTCAT TATCA GAGTA 



CCAGG-3'. The dPiwi flagment was amplified from pGEX-Piwi 
(kindly provided from Dr. Siomi) with the primer set: 5'-GGGCC 
ATGGC TGATG ATCAG GGACG TGGAC GCAGG-3' and 
5'GGGTC TAGAT TATAG ATAAT AAAAC TTCTT TTCG- 
3'. The primer sets were designed based on the registered 
sequence data. The amplified dHPl and dPiwi cDNA were 
digested with jVcoI and Xhol or Xbal and subcloned into between 
the JVcoI site and the Xhol or Xbal site of the pENTRU. The 
resulting plasmids were named dHPl-pENTRl 1 and dPiwi- 
pENTRll, respectively. 
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N-terminal fusion BEST vectors, pIE2-DBD, pIE2-AD, 
pie2RW, pnVW or pcCW, which contains split Venus or 
Cerulean, respectively, were transcribed by the IE2 promoter 
from the Orgyia pseudotsugata nuclear polyhedrosis virus [15,24]. 
The Gateway LR reaction between the entry vectors, BmVLG- 
pENTRU, BmHPla-pENTRll, BmHPlb-pENTRl 1, BmAgoS- 
pENTRll and Siwi-pENTRl 1, and the BEST vectors was 
performed using LR Clonase"™ Enzyme Mix (Invitrogen) 
according to the manufacturer's protocols recommended in 
manufacturer's manual [10,11,25]. The resulting plasmids were 
pBBB-HPla, pBBB-HPlb, pBBB-Ago3, pBBB-Siwi, pBBB- 
dPiwi, pBsRed-VLG, pCC-HPla, pCC-HPlb, pCC-Ago3, pNV- 
Ago3, pNV-Siwi, pAB-HPla and pAB-HPIb. 

The construction of the TetR-based I2H plasmids, the split Piwi 
plasmids for I2H assay and dPiwi expression plasmids is described 
in Methods SI and Table SI. 

Preparation of Double-stranded RNA 

Bouble-stranded RNA (dsRNA) for in vivo RNAi reactions was 
prepared using the following procedures. For the dsRNA template, 
BmAgo3-545 and Siwi-587 fragments were amplified using primer 
pairs: 5'-TCTTC CAAGA GAAGG GTCAA GAAAT AG-3' 
and 5'-CTGTC CCACT AGATA CGAGA GTTTG TG-3'; and 
5'-CCTGA GTTGA TATAT CTAGT GCCAG AAC-3' and 5'- 
TCATA CCTAT CAGCA TAATT CCTAG CC-3', respective- 
ly, from BmN4 cBNA library and were inserted into a blunt- 
ended BsSil-Xbal site of the pENTRl 1. The templates for in vitro 
transcription were synthesized via two-step PGR on the templates 
containing the HPla, HPlb, BmAgo3-545 or Siwi-587 fragments 
using the primers, 5'-TAATA CGACT CACTA TAGGG 
TTTGT ACAAA AAAGC AGGC-3' and 5'-TAATA CGACT 
CACTA TAGGG ACTTT GTACA AGAAA GCTGG-3', for 
the first PGR and the primer, 5'-CAGTG AATTG TAATA 
CGACT CACTA TAGGG-3', for the second PGR. The 
amplified BNA templates amplified were extracted with phenol/ 
chloroform, precipitated with ethanol, and dissolved in H2O. Bi- 
directional transcription reactions were performed in the buffer 
(40 mM Tris-HCl, pH 7.5, 6 mM MgCla, 10 mM BTT, 10 mM 
NaCl, 2 mM spermidine, 2 mM NTPs, 8 units of RNase inhibitor, 
and 20 units of T7 RNA polymerase), and incubated for 2 h at 
37°C. The RNA products were purified and dissolved in 100 mM 
HEPES, pH7.0, incubated for 5 min at 94°C, and left to stand at 
room temperature for 30 min to allow the annealing of the two 
RNA strands, which are used as dsRNA solutions. 

The preparation of dsRNAs against putative piRNA-related 
genes is described in Methods SI and Table SI. 

Cell Culture and Transfection 

The silkworm BmN4 cell line (a gift from Br. Chisa Aoki, 
Kyushu University Graduate School) [11] and BmN4-SIBl 
transgenetic cell lin(; (stored in our laboratory) [17] were 
maintained in IPL-41 medium (Gibco) with W/a fetal bovine 
serum at 27°C. These cells were split at a ratio of 1:7 every 4— 
5 days, seeded in 24-well plates at a density of 5x10 cells/well 
and, one day after seeding, overlaid with a lipid-plasmids complex 
(200 |il/well) in order to undergo transient transfection of 
plasmids. Shortly before the transfection, the lipid-plasmids 
complex was prepared by mixing 1 00 ng of each plasmids with 
5 nl of PBB 1 1 1 solution [26] in a final volume of 30 |4l, followed 
by incubation at room temperature for 20 min and supplemen- 
tation with 170 \A of Sf-900 II SFM (Gibco). The cells were 
incubated for 6 h, maintained for another 72 h in the medium 
replaced with 500 |jl of IPL-41 and the cells were then collected 



for RNA extraction or microscopy. All the images were collected 
using a KEYENCE BZ-8000. 

For RNAi experiments, BmN4-SIBl cells were used, and 
dsRNA for each target gene were added in the medium (100 ng/ 
ml) 3 days before the transfection. To confirm RNAi efficiency for 
each target genes, GAPDH, AG03, SIWI, HPla and HPlb, the 5 
primer pairs: 5'-GGCCG CATTG GCCGT TTGGT GCTCC 
G-3' and 5'-GTGGG GCAAG ACAGT TTGTG GTGCA 
AGAAG-3'; 5'-TCAAT TGACT GGTTT A.\CAG ATGAT 
CAACG-3' and 5'-CAAAC TTTGT TAATG GCCGC GTATC 
TGTCG-3'; 5'-ATGAC TCAGT ACTAT TCAAT CGTCA 
ATGCG-3' and 5'-AACTT GGGAA GAGCA GTTGA AATAC 
AAGTG-3'; 5'-GGTAA AAAAG AGAAG AAAAC GGAGA 
CCAG-3' and 5'-GGCTC GAGTT CTAAT TATTC AGATT 
CGGAT GG-3'; and 5'-GCCGA CAAGA AAAAA GAAAA 
TGAAC CAG-3' and 5'-CCCTC GAGTG TTATT TAGTC 
TTCAT TATTC CC-3', were used respectively. The amplified 
products were separated by electrophoresis through a 1 % agarose 
gel in TAE buffer and stained with ethidium bromide. 

Cell Imaging 

Fluorescence microscopy images were captured using Biozero 
BZ-8000 microscope (KEYENCE). 

Bimolecular Fluorescence Complennentation Assay 

Bimolecular fluorescence complementation (BiFC) analysis was 
based on the reassembling into a functional fluorescent protein 
though the association of protein fragments fused to the proteins of 
interest [13]. 100 ng of each expression plasmid for different 
silkworm Piwi or HPl proteins fused to pnVW or pcCW, 
respectively, were co-transfected into the BmN4 cells. 

Luclferase Assay 

Cells were harvested three days after transfection, washed with 
PBS and lysed with the lysis buffer (25 mM Tris-phosphate, 
pH 7.8, 2 mM BTT, 2 mM Trans-1, 2-diaminocyclohexane- 
N,N,N',N'-tetraacetic acid monohydrate, 10% glycerol, 1% 
TritonX-100). The cell lysate was precipitated at 10,000 rpm for 
1 min at 4°C. Luclferase activity was measured after mixing the 
resulting supernatant and luclferase substrate by ARVO (Perkin 
Elmer) using a luclferase reporter assay system (Promega). 

Supporting Information 

Figure SI The interaction between Piwi and HPl 
proteins in silkworm cultured cells. (A) The interaction 
between D. melanogaster Piwi and HPla in BmN4 cells. The BiFC 
analysis for the interaction between Piwi and HPl proteins in the 
silkworm cultured cell line BmN4. Venus and Hoechst33342 
fluorescences indicate their localization and nucleus, respectively. 
Scale bar: 10 nm. (B) The localization of EGFP-fused Drosophila 
Piwi protein in a BmN4 cell. EGFP and Hoechst33342 
fluorescences indicate their localization and nucleus, respectively. 
Scale bar: 10 \ixa. (C, B) The insect two-hybrid (I2H) assay for the 
interaction between the silkworm Piwi and HPl proteins with the 
use of the TetR-fused Piwi proteins instead of GAL4 BNA- 
binding domain. BmN4 cells were transfectcd with 9xtetO- 
IE2mini(L)-Luc reporter construct, along with ("xpression vectors 
for TetR and p65 activation domain (AB) fused to Ago3, Siwi, 
HPla or HPlb as indicated. Empty cloning vector pUC18 was 
introduced as the control prey instead of pAB-HPla or HPlb. 
Luclferase activities were measured at 72 h post-transfection. 
Values are relative luclferase activities in the pUC18-introduced 
control. Error bars = SB. The SBs and P-values (determined by 
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the Hest, **P<0.01, ***P<0.001, which was checked l)y 
comparing luciferase activities in pUC 1 8-introduced controls with 
pAD-HPla or -b-introduced cells) are based on n = 3. 
(TIF) 

Figure S2 The split Piwi proteins for I2H baits. 

Schematic structure of baits used in the split Piwi I2H assay. Piwi 
protein contains N, PAZ, MID, and PIWI domains. (A) dPiwiFL, 
residues 1-843; dPiwiNT, residues 1-491; and dPiwiCT, residues 
492-843. (B) Ago3FL, residues 1-926; Ago3NT, residues 1-553; 
and Ago3CT, residues 554-926. (C) SiwiFL: 1-899; SiwiNT: 1- 
526; and SiwiCT: 527-899. 
(TIF) 

Figure S3 The N-terminal domain of silkworm Piwi 
proteins interacts more strongly with HPla/b than does 
the C-terminal domain. (A-E) The I2H assay of the 
interaction between the split Piwi and HPl proteins. BmN4 cells 

were transftT,t(;d with 4 xUAS-TATA-Luc reporter construct, 
along with expression vectors for GAL4 DNA-binding domain 
(DBD) and p65 activation domain (AD) fused to dPiwiFL, 
dPiwiNT, dPiwiCT, Ago3FL, AgoSNT, AgoSCT, SiwiFL, 
SiwiNT, SiwiCT, HP la or HP lb as indicated. Luciferase activities 
were measured at 72 h post-transfection. Error bars = SD. The 
SDs and P-values (determined by the t-test, *P<0.1, **P<0.05, 
***P<0.01, which was checked by comparing luciferase activities 
in cells transfected with DBD-fused N-terminal domain plasmids 
with that of DBD-fused C-terminal domain-transfected cells) are 
based on n = 3. 
(TIF) 

Figure S4 The I2H assay for the interaction between the 
silkworm Piwi and HPl proteins in each of silkworm 
Piwi gene knocked-down BmN4-SIDl cells. (A, B) At 72 h 

after the introduction of dsAgo3, dsSiwi or dsGFP (control), 
BmN4-SID cells were transfected with 4xUAS-TATA-Luc 
reporter construct, along with expression vectors for GAL4 
DNA-binding domain (DBD) and p65 activation domain (AD) 
fused to Ago3, Siwi, HPla or -b as indicated. Luciferase activities 
were measured at 72 h post-transfection. The luciferase activities 
were normalized to that of dsGFP-introduced controls (n = 3 
independent luciferase activities). Error bars = SD. P-values are 
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